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1 Introduction 

We are interested in this article in investigating the smoothing effect properties of the 
solutions of the Schrodinger equation. Since the work by Craig, Kapeller and Strausss [13] , 
Kato [IB], Constantin and Saut [12] establishing the smoothing property, many works 
have dealt with the understanding of this effect. In particular the work by Doi [2] and 
Burq [5j [6J shows that it is closely related to the infinite speed of propagation for the 
solutions of Schrodinger equation. Roughly speaking, if one considers a wave packet with 
wave length A, it is known that it propagates with speed A and the wave will stay in any 
bounded domain only for a time of order 1/A. As a consequence, taking the L 2 in time 
norm will lead to an improvement of 1/ A 1 / 2 with respect to taking an L°° norm, leading to 
a gain of 1/2 derivatives. This heuristic argument can be transformed into a proof of the 
smoothing effect either by direct calculations (in the case of the free Schrodinger equation) 
or by means of resolvent estimates (see [3] for the case of a perturbation by a potential 
or [7] for the boundary value problem). In view of this simple heuristics, it is natural to 
ask whether one can refine (and improve) such smoothing type estimates if one considers 
smaller space domains (whose size will shrink as the wave length increases) . A very natural 
context in which one can test this heuristics is the case of the exterior of a convex body (or 
more generally the exterior of several convex bodies), in which case natural candidates for 
the A dependent domains are A _a neighborhoods of the boundary. This is the main aim 
of this paper. To keep the paper at a rather basic technical level, we choose to consider 
only balls, for which direct calculations (with Bessel functions) can be performed. Our 
first result reads as follows: 

Theorem 1.1. Let Q = R 3 \5(0,l), T > 0, < a < § and A > 1. Let and X G C£°(R*) 
be smooth functions with compact support, ijj = 1 near 1, % — 1 near 0. Set X\{\ x \) := 
x(A a (|x| — 1)), where x denotes the variable on Q. Then one has 
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For s G [-1,1] andv{t) = f* e i[ - t ~ T ^ D ip(^)x\gdT 

Xxtp(^f->\\ L ^+\Q) <C^ § \\H^r)X\g\\ L lH h (n), (1-1) 

For s G [0, 1] 



11x^^(^011^.+^ < cxnm^tobw. (i.2) 

ifere i/ie constants C do not depend on T , i.e. the estimates are global in time. 

Using this result, we can deduce new Strichartz type estimates for the solution of the 
linear Schrodinger equation in the exterior, Q, of a smooth bounded obstacle G C M 3 , 

(id t + A D , N )u = 0, it(0, x) = u (x), u\ dn = or d n u \ dn = 0, (1.3) 

that we denote respectively by e !tAD «o and e ltAN uo. 

Definition 1.2. Let q, r > 2, (q,r,d) ^ (2, oo,2), 2 < p < oo. A pair (q,r) is called 
admissible in dimension d if q, r satisfy the scaling admissible condition 

i + ± = i (i.4) 

p q 2 V ; 

Theorem 1.3. (Strichartz estimates) Let 6 = B(0, l)cl 3 and SI = M 3 \ 6, T > and 

(p,q) an admissible pair in dimension 3. Let e > fee an arbitrarily small constant. Then 

there exists a constant C > such that, for all uq G H^ n (Vt) the following holds 

\\e UAD - N u \\ LP{hTyThLqm < C\\u \\ * +£ . (1.5) 

H D,N ( n ) 

Moreover, a similar result holds true for a class of trapping obstacles (Ikawa's example), 
i.e. for the case where is a finite union of balls in R 3 . 

As a consequence of Theorem 11.31 we deduce new global well-posedness results for 
the non-linear Schrodinger equation in the exterior of several convex obstacles, improving 
previous results by Burq, Gerard and Tzvetkov [7]. Consider the nonlinear Schrodinger 
equation on Q subject to Dirichlet boundary condition 

(id t + A D )u = F(u) in 8x0, u(0,x) — u (x), on f2, u\qq = 0. (1.6) 

The nonlinear interaction F is supposed to be of the form F = dV/dz, with F(0) = 0, 
where the "potential" V is real valued and satisfies V(|z|) = V(z), Vz G C. Moreover we 
suppose that V is of class C 3 , \D^V(z)\ < C k (l + |z|) 4 " fc for k G {0,1,2,3} and that 
V(z) > —(1 + l^l)^, for some /? < 2 + |, d = 3 (the last assumption avoid blow-up in the 
focussing case). 
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Some phenomena in physics turn out to be modeled by exterior problems and one may 
expect rich dynamics under various boundary conditions. A first step in this direction is 
to establish well defined dynamics in the natural spaces determined by the conservation 
lows associated to (11.61) . If u(t, .) G Hq(Q) f]H 2 (Q) is a solution of (11. 6ft then it satisfies 
the conservation lows 

f \u(t,x)\ 2 dx = 0; [ \Vu(t,x)\ 2 dx + [ V(u(t, x))dx) = (1.7) 

dt Jq dt \ Jq J n J 

and therefore for a large class of potentials V the quantity \\u(t, -) || j?i(n) remains finite along 
the trajectory starting from uq G Hq(Q) f]H 2 (Q). This fact makes the study of (II. 6p in 
the energy space Hq(Q) of particular interest. It is also of interest to study (11.61) in L 2 {Vt): 
the main issue in the analysis is that the regularities of H 1 and L 2 are a priori too poor 
to be achieved by the classical methods for establishing local existence and uniqueness for 
( II. 6p . We state the result concerning finite energy solutions, which will be a consequence 
of Theorem 11.31 

Theorem 1.4. (Global existence theorem) Let Q = M 3 \ 5(0,1). For any uq G Hq(Q) 
the initial boundary value problem (11. 6ft has a unique global solution u G C(M., Hq(Q)) 
satisfying the conservation laws (11.71) . Moreover, for any T > the flow map uq — > u is 
Lipschitz continuous from any bounded set of Hq{VL) to C(R, Hq(Q,)). 

Remark 1.5. In [7], Strichartz type estimates with loss of - derivative have been obtained 
for the Schrodinger equation in the exterior of a non-trapping obstacle C M. d which 
allowed to prove the same global existence result in dimension 3 provided ||wo||jfi(n) is 
sufficiently small. 

The Cauchy problem associated to (11.61) has been extensively studied in the case Q = 
M. d , d > 2, by Bourgain [I], Cazenave [10], Sulem et Sulem [23J, Ginibre and Velo [15], 
Kato [16] and the theory of existence of finite energy solutions to (11. 6p for potentials V 
with polynomial growth has been much developed. Roughly speaking the argument for 
establishing finite energy solutions of (11. 6p consists in combining H l local well-posedness 
with conservation laws (I1.7P which provide a control on the H 1 norm. 

The article is written as follows: in Section [2] we obtain bounds for the L 2 norms for 
the outgoing solution of the Helmholtz equation that will be used in Section [3] in order to 
prove Theorem 11.11 In Section H] we use a strategy inspired from [22] , [7] to handle the 
case when the solution is supported outside a neighborhood of dQ; in Section [5] we achieve 
the proof of Theorem 11.31 The last section is dedicated to the applications of these results; 
precisely, we give the proofs of Theorems 11.31 and 1 1 . 41 in the case where the obstacle consists 
of a union of balls. In the Appendix we recall some properties of the Hankel functions. 

Acknowledgments: This result is part of the author's PhD thesis in preparation at 
University Paris Sud, Orsay, under Nicolas Burq's direction. 
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2 Precise smoothing effect 
2.1 Preliminaries 

Let Q C M d , d > 2 be a smooth domain . For s > 0, p G [1, oo] we denote by W' p (f2) 
the Sobolev spaces on Q. We write L p and .fP instead of W 0,p and jy s ' 2 . By Ad (resp. 
A = An) we denote the Dirichlet Laplacian (resp. Neumann Laplacian) on Q, with domain 
H 2 (n) n H^(Q) (resp. {u G # 2 (fi)|<9„u| m = 0}). We next define the dual space of i^(ft), 
iJ _1 (f2), which is a subspace of D'(Q). We construct H~ S (Q) via interpolation and due to 
[21 Cor.4.5.2] we have the duality between H$(tt) and H~ s (£l) for s G [0, 1]. 

Proposition 2.1. Let d > 2 and Q C M> d be a smooth domain. The following 
continuous embeddings hold: 

• /To 1 (ft) C L 9 (fi), 2<q<^ (p < oo if d = 2), 
. C 1-1 = ^, aG [0,1), 

. HFMCW^M, i-l = £, a 6[0,l), 

• W^(fi) C L°°(ft) ; s > |, p > 1. 

Proof. The proof follows from the Sobolev embeddings on M d and the use of extension 
operators. □ 

Let ip,x G C^°(R*) be smooth functions such that x = 1 in a neighborhood of and 
for all r G R+, E fc >o^( 2 ~ 2fcr ) = L For A > let Xa(M) := x(A Q (M - 1)), where a; 
denotes variable on Q. We introduce spectral localizations which commute with the linear 
evolution. Since the spectrum of -A^ is confined to the positive real axis it is convenient 
to introduce A 2 as a spectral parameter. We will consider the linear problem (jl.3p with 
initial data localized at frequency A, u|t=o = ^r) u o- 

Remark 2.2. In order to prove Theorem [L3] it will be enough to prove (jl.5p with the initial 
data of the form tf}(— jp)uo, since than we have for some e arbitrarily small 

||e ftAD « |k P ([-T,T],w(n)) ~ ll(e i * AD ^(-^)«o)il|w([-T,T] > i«(n))z| < (2- 1 ) 



M 4+ £ • 

^ («) 
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2.2 Estimates in a small neighborhood of the boundary 

In what follows let d = 3. We study first the outgoing solution to the equation 

(A D + X 2 )w = X xf, w\ dn = 0. (2.2) 

In what follows we will establish high frequencies bounds for the L 2 norm of w in a small 
neighborhood Q\ = {\x\ < A - "} of size X~ a of the boundary. We notice that a ray with 
transversal, equal-angle reflection spends in the neighborhood Q\ a time ~ A~ Q . If the ray 
is diffractive then the time spent in Q\ equals A". 

We analyze the outgoing solution of (12. 2p outside the unit ball of M 3 . The first step 
is to introduce polar coordinates and to write the expansion in spherical harmonics of the 
solution to 

(A D + X 2 )w = on tt = {x G M 3 ||a;| > 1}, 
15180 = / on § 2 , (2.3) 
r(d r w — iXw) — > r _»oo 0. 

In this coordinates the Laplace operator on Q writes 



A D = d 2 + -d r + \A § 2, (2.4) 



where A S 2 is the Laplace operator on the sphere S 2 . Thus the solution w of (12.31) satisfies 

r 2 d 2 w + 2rd r w + {X 2 r 2 + A S 2)w = 0, r > 1. (2.5) 

In particular, if {ej} is an orthonormal basis of L 2 (S> 2 ) consisting of eigenfunctions of A§2, 
with eigenvalues — /i 2 and if u denotes the variable on the sphere § 2 , we can write 



wiruj 



J2^(r) ej (u), r>l, (2.6) 



where the functions Wj(r) satisfy 



r 2 w'^(r) + 2r^.(r) + (A 2 r 2 - f£j)wj(r) = 0, r > 1. (2.7) 

This is a modified Bessel equation, and the solution satisfying the radiation condition 
r(d r w — iXw) ^ r ->oo is of the form 



w j (r)=a j r-^H VJ (Xr), (2.8) 
where H v . (z) denote the Hankel function. Recall that the Hankel function is given by 

9 i(z-7vu/2-ir/4) roo 

^ = (^ 1/2 Vrwi '-•r-^l-^r**. (2.9) 

where 

/■oo 

T(z/+l/2)= / e~ s s u - ll2 ds, 
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and it is valid for Rev > | and —tt/2 < argz < ir. Also, in (12.81) Uj is given by 

"j = (ti + \) 1/2 (2-10) 
and the coefficients a,j are determined by the boundary condition Wj(l) =< f,ej >, so 

a = < J' e /> . (2.11) 



Let us introduce the self-adjoint operator 



A = (-Aga + ^) 1/2 , Ae,- = z/j-e,-. (2.12) 



Then the solution of (12.31) writes, formally, 

^M=r-^^M/ M , WG§2 . (2.13) 

Proposition 2.3. (see Chp.3]) The spectrum of A is spec(A) = {m + -Am G N}. 

Remark 2.4. (see [211 Chp.3]) The Hankel function H m+ i/ 2 (z) and Bessel functions of order 
77T.+ I are all elementary functions of 2. We have 



H m+1/2 (z) = (-) 1/2 q m (z), q m (z) = -i(-l)™(ii!-)™(£_). (2 .14) 

7T 2 z 



We deduce that 



r _ 1/2 tf^_) = r . m . Ie , A ,. 1) p^_) i = ._„_, £ + 

F m +i/2(A) Pm(A) ^ fc!(m-fc)! 

(2.15) 

We now look for a solution to ( 12.21) . If we write 

f(ru) = /i( r )eiH, w(ru) = ^ Wj(r)ej(u), 
j j 

then the functions Wj(r) satisfy 

r 2 w"(r) + 2rw' j (r) + (A 2 r 2 - fi 2 )wj{r) = r 2 fj(r), r > 1 (2.16) 

together with the vanishing condition at r = 1 and Sommerfield radiation condition when 
r — > oo. Applying the variation of constants method together with the outgoing assumption 
and the formula (12 . 1 3j) . we obtain 

roc 1 

w,(r)= G I/j (r,s,X)xx(s)f 3 (s)s 2 ds, v s = (tf + -) 1/2 , (2.17) 
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where G u (r, s, A) is the Green kernel for the differential operator 

^ = -^ + --f + (A 2 -^)> -=(^) 1/2 . (2.18) 
ar r ar r z 4 

The fact that L v is self-adjoint implies that G u (r, s, A) = G u (s, r, A) and from the radiation 
condition we obtain (for A real) 

|(r S )-^(j„( r A) - £$/r„(rA))j/„(«A), r < >. 

In what follows we will look for estimates of the L 2 norm of Wj on the interval [1,1 + A~ a ]. 
This problem has to be divided in several classes, according whether Vj/X is less than, 
nearly equal or greater than 1. We distinguish also the simple case of determining bounds 
when the argument is much larger than the order. Let us explain the meaning of this: in 
fact, applying the operator L v . to rwj(r) instead of Wj(r), we eliminate the term involving 
w'j(r) in (12.161) . On the characteristic set we have 

P] + ^ = W ^ = (/.? + i) 1/2 , (2.20) 

where pj denotes the dual variable of r. Let 9j be defined by tanfy = then for A big 
enough and r in a small neighborhood of 1 we can estimate 

A 2 

tan 2 9,-, + 1 ~ 

• When the quotient y is smaller than a constant 1 — e where eo is fixed, strictly 
positive, this corresponds to an angle 9j between some fixed direction 9q and tt/2, 
with tan#o = e o ; and thus to a ray hitting the obstacle transversally. In this case we 
show that, since a unit speed bicharacteristic spends in the A-depending neighborhood 
Q\ a time X~ a , we have the following 

Proposition 2.5. Let eo > be fixed, small. Then there exists a constant C = C(eo) 
such that 

ll w illL 2 ([l,l+A- 0! ]) < CX (1+a) ||/j||L 2 ([l,l+A- Q ]) 

uniformly for j such that {vj/\ < 1 — eo}. 

• When the quotient y is close enough to 1 the angles 8j become very small and this 
is the case of a diffractive ray, which spends in £l\ a time proportional to A~ Q / 2 . In 



this case tan 9j ~ y 1 — and we show the following 
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Proposition 2.6. If 1 — -f ~ A ^ /or some /? > £/ien 

INIr»([i,i+A-]) < ^-^"^^^^([Li+A-]) z/ /3 < a 

and 



||^IU- C [i,i+a-«D < CX~ {1+a/2) \\f 3 \\ LH[1 , 1+x -^ ifP> 



a. 



• In the elliptic case ^- C 1 or ^- 6 [eo, 1 — eo] f° r some small eo > there is nothing 
to do since away from the characteristic variety (I2.20p we have nice bounds of the 
solution of the solution Wj(r) of (12.161) . 

Proof. ( of Proposition [O]) The solution Wj(r) given in (12 . 1 7[) writes 

H Vj {\) 



(r) = J\x{s)ms" 2 {H Vj {\s) - J ^H Vj (Xs))dsH Uj (Xr)+ [2.21] 

/DO g ( \ \ > 

Xx {s)f 3 {s)s^H U] {\s)ds{H Uj {\r) - -^. Hv ,{\r)) ^ 



thus in order to obtain estimates for (t*) ||x^([i,i+a— «]) we have to determine bounds for 
||/J^(Ar)|| i 2([ 1 1+ A-aj) since we have for some constant C > 

lkiWIU 2 ([i,i+A-]) < \\fj(r)\\ L 2 {[1:1+x - a]) \\H U3 (Xr)\\l 2m+x - a]) . (2.22) 

We consider separately two regimes: 

1. For !f < 1 we use (17. 2p to obtain immediately 

ll^i(Ar)||i. aifl+A - ]) <A- 1 - a . (2.23) 

2. For < e < Cj := y < 1 — e for some fixed e > we use the expansions (17.31) . 
(Haj): set 

Cj sin r c 7 - 

cos r = — , dr = c,- 



— — rfr, r G [arccoscj, arccos - — ^_ a ] ='■ [ 7 j,o>'0',i] 
Then we have 



3 cos^ 



J u ' J; ttA /_.. cos t 7TA Vl-sinr/ J: 

1/ (1 + Jl-^/(l + A-)2) 2 ! A -a 

~ —[ (1+A^ a ) 2 " -== 1 ] ~ - - 

" TTA V 1 ^ (1 + ^13^2)2 / ~ TTA (1 - C 2 + A"") 1 / 2 + (1 - C 2 )V2 

and 1 — c 2 G [2e — e 2 ,, 1 — e(j] with < e < 1 fixed. Notice, however, that if one 
takes Cj ~ A _/3 for some /? > 0, then tanflj ~ (1 — c 2 ) 1 / 2 ~ A _/3 / 2 . If /? < a we can 
estimate f[2~24l) by A-CWa-^-/^ otherwise we have the bound A" (1+a/2) . 
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□ 

Proof, (of Proposition \2.6§ Here we use Proposition 17.31 Let 1 — Uj/X = tX^ 13 for r in a 
neighborhood of 1. For s e [1,1 + X~ a ] write z(s)uj = sX, thus z(s) = 1 _ 7 A A _ )3 . Since 

A -a/2 

IKIUa([i,i+A-a]) < C ^ H ([1,1+a-"]) 1 1 H Uj ( Ar ) 1 1 L 2 ([i,i+A~ a ]) x (2.25) 

x||J i/j (As)n j (A) - J, 3 (A)n j (A S )|| L2(1 , 1+A - a) , 
we shall compute separately each factor in (I2.25P (modulo small terms). We have 

ll^(Ar)||| 2(1)1+A - a) = |^(Ar)| 2 + \Y Vj (Xr)\ 2 dr ~ (2.26) 

2 _j /- 1+A ~ a 1 2 _ x X~ a / A- 1 -^-^), si 



tt^' A (^ 2 (s) - iy/ 2 ds ~ ir Uj X- a + 2tA-^ ~ \ A" 1 , si, p >a, 

\H„M)\ 2 = \^W\ 2 + I^WI 2 - ^7 1 ( ^ 2(1) 1 _ 1)1/2 - A~ 1+/3/2 , (2.27) 
while for the factor in the second line in (I2.25P we have 

\\J Uj (Xs)Y Uj (X) - J Vj (X)Y Vj (Xs)\\ L 2 (1>1+x - a) < (2.28) 



n< {z 2 {l)-l) 1 /^J 1 (z 2 {s)-iy/ 2 J -\A- 1+ ^ 4 , si, (3>a. 
From fl2T25|) . (I2T26D . fl2T27]) . (I2T28D we deduce 



A -l-a/2-(a-/3) j g j < /3 < tt, 
j||L2([l,l+A— ]) < C||/j||L2([l,l+A— ]) X <j A _l- a /2 gi ' p > a ~ ( 2 ' 29 ) 



□ 

Neumann: As far as the Neumann problem is concerned, we must solve the problem 



^ 2 + + (A 2 - ^j)v>j(r) = x(A°(r - l))/,(r), -^(1) = 0, (2.30) 
which gives, after performing similar computations 



- k~ 1/2H ^ Xr K]HfW 2 I ^ s ^> 3/2 H^)ds- (2-31) 



r 



poo 

Propositions \2. 5[ \2.6i hold true for the Neumann case too. 



3 PROOF OF THEOREM 1 



10 



3 Proof of Theorem 1 

3.1 Reduction to the Helmholtz equation 

Proposition 3.1. Consider the Helmholtz equation 

(A D + X 2 )w = X xf on n 



(3.1) 

w\an = 0, 

with the Sommerfield "radiation condition" (where r = \x\) 

r(d r w — i\w) — > r ^oo 0. (3.2) 

In order to prove Theorem \1.1\ it is enough to establish estimates for the L 2 norms of the 
Helmholtz equation (13.11) . 

Proof. Consider the inhomogeneous Schrodinger equation 

(id t + A D )v(t,x) = x\9(t,x), v\ t=0 = 0. (3.3) 

We denote by v ± the solutions to the equations 

(id t + A )^(f,i) =x\g(t,x)l { ± t>Q y, v ± \ t= o = 0. (3.4) 

For e > and ±t > we define vf = e^ et v which satisfy 

(id t + A D ± ie)vf = l{± t >o}e~ et X\g(t, x), ^ ^ 
vt\t=o = 0, |an = 0. 

After performing the Fourier transform T with respect to the time variable t, the equation 
(13.51) becomes 

(A D - t ± ie)vf(r, x) = X\^(^-{±t>o}e^ et X\g)(T, x), 



f;± , _ n (3-6) 
v e |an — U, 

where r denotes the dual variable of t and we deduce 

Xa£+(t, x) = X\(A D -r + iey 1 X\J r (l{±t>o}e~ et X\g)(T, x). (3.7) 

Since -A^ is a positive self-adjoint operator, the resolvent (—Ad — z)" 1 is analytic in 
C\ R + . Since the spectrum of — is confined to the positive real axis it is convenient to 
introduce A 2 G M + as a spectral parameter. Notice, however, that there are two manners 
to approach A 2 > in C \ M. + , choosing the positive imaginary part, which corresponds 
to considering A 2 + ie, or the negative imaginary part, which corresponds to A 2 — ie. The 
"physical" choice corresponds to the limiting absorption principal and consists of taking 
A 2 + ie. In some sense, the limiting absorption principal allows to recover the "sense of 
time" . If one replaces A 2 by A 2 ± ie, e > 0, then (A 2 ± ie) belongs to the resolvent of the 
Laplace operator on Q with Dirichlet boundary conditions on dQ and we can let e tend to 
in (13.71) (since the operator xx(Ad — t + ie)~ l x\ has a limit as e — > 0) and so we can 
express the Fourier transform of the unique solutions v ± of (13. 4p as 

Xa^ ± (t,x) = \imxx(A D - t ±iey 1 x\J r (l{±t>Q}^ tt Xxg)(r,x). (3.8) 
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Remark 3.2. Notice that here we used the fact that for — p- away from a neighborhood of 
1 we have the bounds 

ii H-H ^ C\\x\g\\ L 2 
Wxxv \\ L2 < -y^- 

We conclude using that if w e = t>+, f e = !F{l{t>o}er €t X\9) the following holds 

Proposition 3.3. (124-1) As e — > 0, w e converges to the unique solution of f13.in - fj3.2p . 
where f = lim e ^ fe- 
lt remains to notice that if 7i is any Hilbert space, than the Fourier transform defines 
an isometry on L 2 (IR, 7i). □ 

3.2 Smoothing effect 

In this section we prove Theorem 11.11 

• We prove ( 11 .11) for s — 0. Let A > 0, w and g be such that (12.21) holds. We multiply 
( 12.21) by xxw and we integrate on fl 



Xx\Vw\ 2 dx = X 2 j x\\w\ 2 dx-2\ a J < Vw, Vx(A a (|a:| - 1)) > wdx- J xlfwdx, 

(3.9) 

thus, using that < xa < 1, Xa < Xa> Vx(A a (|x| — 1)) < Xa, we have for all 5 > 

J \x\^w\ 2 dx < X 2 J \x\w\ 2 dx + 5 J |xa/| 2 ^+^ J \x\w\ 2 dx + \ a j w\Vw\x\dx. 

(3.10) 

Since for all 8\ > one has 

J w\Vw\ X \dx < ( J \ X xVw\ 2 dx) 1/2 { J \ X xw\ 2 dx) 1/2 < ^XaVHII^+^IIxaHII^), 

taking 5 = A~ Q , 6\ = \~ a /2 together with J |xa^| 2 ^ ^ 2a I \Xxf\ 2 (according to 
the computations made in the preceding section) we deduce 

J \ X xVw\ 2 dx < \~ a J \ X xf\ 2 dx. (3.11) 

Thus we have obtained 

IIxaHI^) % A~f ||xa/||l2(o). (3.12) 

Dualizing (13.121) we find 

IIxaHU 2 ^) < A"f \\x\f\\ H -\ny ( 3 - 13 ) 
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which will yield (11.11) for s = — 1. Now we prove (11.11) for s = 1. Let again w and 
/ be such that (12. 2p holds and let xx be a smooth cutoff function equal to 1 on the 
support of xx- Write 

\\Xxw\\ H 2 D{n) « || x\w | \ H i^o) + ||A D (xAw)|| L2(n) . (3.14) 

Since HxA^H^tn) can be estimated by means of (I3.12p . we only need to obtain 
bounds for || A£>(xaw)||l 2 (c)- We write 

A D (xxw) = XxA D w + [A D , Xxlxxw. (3.15) 

The commutator [Ad, xx] is bounded from H D (Ji) to L 2 (Q) with norm less then C\ a 
and we have chosen a < 1 (C does not depend on A), while the first term in the 
right hand side of (13.151) is in H D (Q) since Apto = xxf — A 2 w and satisfies (12. 2p with 
Xxf replaced by Ad(xa/)- Therefore, we can apply (13.121) in order to deduce the 
inequality (11.11) for the Fourier transforms in time of u and / which appear in the 
proposition. Since we have obtained the result for s = — 1 and s = 1 we can use an 
interpolation argument to get it for s G [—1, 1]. 

• We turn to the proof of (11.21) for s = 0. If we denote by A\ the operator which 

to a given u e L 2 (f2) associates Xxe ltAD ip(-j^-)uo, we need to prove that A\ is 

i 

bounded from L 2 {Vt) to L^H^iVL) with the norm less than C\~^ for some constant 
C independent of A, which in turn is equivalent to the continuity of the adjoint 
operator, 

AUf) = j%(^)e- A -Xxf(r)dr, (3.16) 

from L\H D 2 (VL) to L 2 (f2) with the norm bounded by C\~^, which is equivalent to 
showing that the operator A X A* X defined by 

(A x Alf)(t) = j\ xe ^^(^)e- A -xxf(r)dr (3.17) 
_ i i 

is continuous from L\H D 2 (VL) to L^H D (Q) and its norm is bounded from above by 
CA-2. We write (A x A\f)(t) as a sum 

(A x A*J)(t) = j\^{Z^yit-r)^ { Z^ )xJ{j ) dT+ ( 3.i 8) 

+ f Xx^^)e^ A ^(^) X xf(r)dr. 

Hence, in order to conclude it is sufficient to apply (11.11) with s — — ~ together with 
time inversion, since the second term on the right hand side of (13.181) will solve the 
same problem with initial data u\ t =T = uq- 

We prove now (11.21) for s = 1. 
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Lemma 3.4. We have 

\\ A ^ eltAD ^^>o)\\ L2TH -j {n) < X^WxxH^HWmny (3.19) 



Corollary 3.5. Lemma 3Jf_ yields 

IIXAe l4A ^(^)n || L ^| (n) < X^WxM^UoWl^). (3.20) 

Corollary 13.51 and an interpolation argument now yield for 9 G [0, 1] 

IIXAe l4A ^(^) Mo |l L ^| (c) < A"-^||xa^(-^KIU 2 (q), (3.21) 
achieving the proof of Theorem 11.11 
Proof. ( of Lemma \3.J$ Write 

(-A D + 1) ( XA e^(^) Mo ) = Xx (- A D + 1) (>^(^H) - (3.22) 

-[A D)XA ](e^^(^) Mo ). 
— For the first term in the right hand side of (13.221) we show that the operator 

B x := ^(fi)9 % ^ XA (-A D + l)( e itA ^(^H) e4ff^(n)) (3.23) 

_i 

is continuous and its norm from H})(Q) to L^H D 2 (Q) is bounded from above 
by CA - ^ for some constant C independent of A, or equivalently that 

[b x (-A d + l)- 1 ^/)^) = J\x{-^d + l)^ 2 ( Z ^)e i( '- T)A -XA/(r)rfr = 

= (-A D + l)(A A A*/)(t) + [A D , Xx ] J T e^ ) A D ^ {Z ^ )xxf{r)dr 

(3.24) 

is bounded from L^H^(Q) to L\H D 2 (Q) by CA _ 2\ Here is the operator 

introduced in the proof of the case s = 0. For the first term in the right hand 

side of (13.241) we apply (11.11) with s — ~ and we obtain a bound CA~ 2, while for 

the second term we make use of (II. 2p with s = — \ and of the fact that [An, x\] 
1 _ 1 

is bounded from H D (Q) to H D 2 (Q) with a norm bounded by CA a in order to 
find a bound from L\H D {Vt) to L\H D 2 {Vt) of at most CA^ -1 < C\~%. 
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The second term in the right hand side of (I3.22p gives 



"An . ,, _.„,,,, / _ , , — A; 



For evaluating the first term in the last sum we use again (11.21) with s = 

A-Ke«^(^) Mo || , < A-f ||^(^£) Mo | Ua(!1 , (3.26) 



where xa is a smooth cutoff function such that xx is equal to 1 on the support 
of xx an d we conclude since a < 1. For the second term we have 



< A 2 | x'a 2 | (e^(=^)«o)«| + 2A« +1 / 2 ||xA||i 2( n) < A 2 ||xa^||| 



(3.27) 



where we have set "0(x) = xip(x). From (11.21) with s = that we have already 
established, we find, since a < 1, 

mv^^H)!!^^ < A1Ml i|Hj(n , (3.28) 

^ ,/~ A/? x i_o M „ Ac. 

^ A * \\x\V{—jp-)Uo\\v(n) ~ A 4 IIxa^(^-)mo||l2(c)- 

□ 

Neumann: v4Z/ of the above results remain valid if we consider the Neumann Laplacian 
An: in fact let w and / be such that (12.21) holds. Using the same strategy as before we get: 



d n u(xxu)da— xa|Vw| dx+X / xx\ w \ dx—X a / < Vw, V%a > Xxwdx = / Xxf w dx- 
on Jn Jn Jn Jn 

(3.29) 

Notice that, due to the Neumann boundary conditions, the first term in the left hand side 
vanishes, so the computations are almost the same as in the previous case. However we 
must pay a little attention to the spectrum of the Neumann Laplacian: for, we have to 
introduce a spectral cut-off <p G equal to 1 close to and decompose 

w = (f>{-A N )w + (1 - <f>{-A N ))w, f = <P(-A N )f+(l-<p(-A N ))f, (3.30) 

m = (p(-A N )u + (1 - <p(-A N ))u . 

We can then rewrite the above proof with w, f , Uq replaced by (1 — </>(— An))w, (1 — 
A^))f, (1 — <J)(—An))u in order to obtain similar estimates in Theorem \l.l\ To deal 
with the contributions of the remaining terms we use the fact that in this situation the L 2 
and H N norms are equivalent. 
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4 Estimates away from the obstacle 

In this section we obtain bounds away from the obstacle. The main idea is to construct 
a new function v = <pu which will solve a problem with a nonlinearity supported in a 
compact set; under these assumptions, it is proved that the free evolution satisfies the 
usual Strichartz bounds (see [22] ). However, we have to take into account the fact that the 
neighborhood outside of which we will apply this result is of size X~ a and thus we will "lose" 
a derivative; however, this will not pose any problem if a is chosen small enough, since it 
will be covered by the loss of derivatives near the boundary. After a change of variables 
we can assume that Q = {x 3 > 0} and thus x 3 defines the distance to the boundary. We 
define e C°°(0) by 0(aOU s <i = x 3 , (f)(x)\ X3 > 2 = 1. 
Proposition 4.1. We have 



1 - Xx)e ltA ^(^)u \\ L ^ (n) < X a \m^)u \\mn } . (4.1) 



We set v = (jm, where u = e ltAD ip(—£f-)uo solves (jl.3p . Then v solves the equation 



f (idt + A D )v = (A D( j))u + 2V0Vu + g^ n ® - v\ dn ® tfgj = [A D , </>]u, ^ 2) 
\ v\ t=0 = (f)u\ t=0 , v\ dn = 0, 

where 6 is the Dirac measure on dQ. It can easily be seen that the last two terms vanish. 
We have [A D ,(p]u G L 2 T Hj mp (VL) and (see Prop.2.7]) 

(4.3) 

The inhomogeneous part of the equation (14. 2 p satisfied by f has compact spatial support 
and therefore we can employ [221 Thm.3] (for p = 2) and Prop. 2. 10] in order to obtain 
Strichartz estimates without losses for v 

ll0e l4A M^HIU^(c) < II^^HIU^), (4.4) 
for every (p, q) <i-admissible pair, which in turn yields 

A— - XaK^^NIIl^) < 11(1 - XxW tA ^(^)uo\y TL ^) (4.5) 



In particular, for p = 2 we get (14.11) from (14. 5 j) . 

Neumann: When dealing with the Neumann problem, this approach requires some 
adjustments, but the above result remains valid (with a slightly modified). Let e > and 
consider <f>\ xa <i{x) = x\ +e , <^|a;3>2 (^) = 1- Then v = (pu solves the equation 

(id t + A D )v = [A D ,(p 1+e ]u, v\ t =o = <j>u\ t =o, -^\an = (x € 3 u)\ 9n + (p—\ 9n = 0. (4.6) 

It's easy to see that one can obtain similar estimates, with a replaced by a — e; still, this 
makes no difference for our purpose, since e can be chosen as small as we like. 
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5 Proof of Theorem 2 

In this section we achieve the proof of Theorem 11.31 Taking s — 1/2 in (11.2j) gives 

\\ Xx e itA °^(^^)u \\ L ^ Hhm < A^-f ||xa'0(^^)«o|U 2 (Q) (5-1) 



from which we deduce 

\\x^ tA ^(^)uo\\ L2L ^ m < \\xxe^(^)u \\ LlHh(n) < X^U(^)u \\ LHU) . 

(5.2) 

An energy argument yields 

\\X\e lt ^ D ij{-^f-)uQ\\ L ^ L 2 {n) < ||^(_^£)u ||ia(n). (5.3) 
Interpolation between (15. ip and (15.31) with weights - and 1 respectively yields 

\\Xxe UA ^(^)u \\ LPTLHn) < A^ 1 -t)||^(^HHz, 2 (Q). (5.4) 
We have also obtained estimates away from the boundary 

||(1 - XA y* A ^^HIU* W (Q) < \ a M^)u \\ L i [n) . (5.5) 

If we take a = -(1 — |) we find a = Now, for p = 2 this gives a = |(< |) and 

consequently we have 

l|e" A ^(^HII L ^ (c) < Xhm^Hhnny (5.6) 
Interpolation between (I5.6P and 

]le uA D ^ {Z ^E )uoh¥L2{n) < ||^ ( Z^ )uo || L2(n); (5.7) 

with weights ~ and 1 — 2 yields, for every (p, g) admissible pair, p> 2 



|e ftAD V(-^ £ )«o||i5.L.(n ) < A*IIiK-^)«oIIl»(0). 



Neumann: The case of the Neumann conditions could be handled in the same way. 
However in (15.51) we have a — e instead of a so we find a = | + and for every (p, q) 
c/-admissible pair and every e > we have 



\\j tAN i>(^)u \\L» T L, m < \^m^)uo\\mn). (5.9) 
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6 Applications 

6.1 Proof of Theorem 11.41 for non-trapping obstacle 

The proof of Theorem 11.41 relies on the contraction principle applied to the equivalent 
integral equation associated to (11.61) with Dirichlet boundary conditions 

u(t) = e itAD u - i [ e i{t - T)AD F(u(T))dT. (6.1) 
Jo 

The assumptions on F and on the potential V imply the following pointwise estimates: 

\F(u)\ < \u\(l + \u\ 2 ), \VF(u)\ < \Vu\(l + \u\ 2 ), \F(u)-F(v)\ < \u-v\(l + \u\ 2 + \v\ 2 ), 

(6.2) 

\V(F(u) - F(v))\ < |V(u -v)\(l+ \u\ 2 + \v\ 2 ) + \u- v\(\Vu\ + |Vv|)(l + \u\ + \v\). (6.3) 

The aim is to show that for sufficiently small T > we can solve the integral equation 
by a Picard iteration scheme in the Banach space Xt ■— L^Hq(Q) fl L^W cr,q (Q), where 
2 < p < y and a — 1 — ^. We equip X T with the norm 

IMUt := IMU°?flJ(n) + 11(1 - A D )?n|| L P L9(c) . 

Remark 6.1. For 2 < p < ^ one has the continuous embedding W l ~^p ,q {Vt) C L°°(f2) 
where (p, q) is an admissible pair, since in this case 1 — ^ > | and Proposition 12.11 implies 
u E L T L°°({1) for any u £ Xt. 

Define a nonlinear map $ as follows 

$(«)(*) : = e^Mo - 2 / e i( '- r)AD F( M (r))rfr. (6.4) 

We have to show that $ is a contraction in a suitable ball B(0,R) of Xt. Using the 
Minkowski integral inequality together with an energy argument we have, on the one hand 

ll^Wllx-ffi(n) % ll«o||flj(n) + H^WIU^n) ^ 

, . || 1 1 1 1 1 1 1 1 2 / II ii / 1 ~ 1 1 1 1 1 1 1 1 2 \ 

(pollen) + ll M llL^ ( i(n)l|M|| L? >^i (Q) + T\\u\\ L¥H i {Q) +T v \\u\\ L¥H i^\\u\\ l p loo( ^) < 

c(\\u \\ HU n) + T\\u\\ Xt + (T 1 -! + T)\\u\\\ T ) (6.5) 

and also 

||(l-A z ,)f<5( w )|U^ (n) < (||(l-A z ,)f^ || ^_+/ T ||(l-A D )f(l+Ks)| 2 M S )|| * )< 
(IMIfl3(n)+jf (l+\Hs)\\l~)\Hs)\\ H i {Q) ) < c(||wo|| if i (n) +r||u|| XT +T 1 ~i||u||^. (6.6) 
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We chose R > AC\\u \\ H i^ and T such that AC\T + (T + T p)R 2 j < 1. In a similar way 
we obtain 

||$(u) - || Xt < C||« - w||x T (T + (T + T 1 -I)(l + ||«||^, + \\v\\ 2 Xt )) . (6.7) 

Taking T eventually smaller such that c(t + (T + T 1 ~p)(l + 2R 2 )j < 1, we deduce from 

fl63|) . m and (Q that $ in a contraction from 5(0, 5) C X T to 5(0, R). Therefore, if 
we consider the sequence {f n }neN C Xt such that v o = uq G B(0,R), v n+ i = &(v n ), then 
f n converges in Xr to the unique solution in Xt of the integral equation 

u(t) = e ltAD u - i [ e i{t - T)AD F{u{r))dr (6.8) 
Jo 

which yields the local well-posedness result. From (16.71) we obtain the Lipschitz property 
of the flow map. Using a standard approximation argument we can derive the conservation 
laws. Next due to fll.7p . the assumption on V (V(\u\ 2 ) > 0) and the Gagliardo-Nirenberg 
inequality we can extend the local solution to an arbitrary time interval by reiterating the 
local-posedness argument. 

Neumann: When we consider the Neumann Laplacian we take Xt like before and 
we choose e so small such that the embedding W 5 p 5 p' q (tt) C L°°(f2) still holds: for 
e < 12 ~ 5p we are led to the same conclusion. 

lb 

6.2 Proof Theorem 11.31 for a class of trapping obstacles 

In this part we prove Theorem 11.31 for a class of trapping obstacles. In this case, since there 
are trapped trajectories (e.g. any line minimizing the distance between two obstacles has 
an unbounded sojourn time), the plain smoothing effect Hz does not hold. However, one 
can obtain a smoothing effect with a logarithmic loss (see [H Thml.7,Thm.4.2]). 

Assumptions: We suppose here O = U^ 1 0j is the disjoint union of a finite number of 
balls 0j = Bi(oi, Ti) of radius rj > in M 3 . We denote by k the minimum of the curvatures 
of the spheres §(r^) = <90j, that is k = mm{^=,i = 1,N}. Denote by L the minimum of 
the distances between two balls. Then, if N > 2, we assume that kL > N (L > I > if 
iV = 2 for some strictly positive /). We keep the notations of the previous sections. 

Let x ^ 1)), x — 1 close to and for every i, set Xi( x ) '■= x( ^7°^ — -0 anc ^ 

Xi,\( x ) — x(^ a (^7^ ~~ 1)) (Xi,x vanishes outside a neighborhood of size \~ a of the ball 
0j). Set u = e ltAD i/j( :L ^-)uo and Vi = Xi u - Then Vi satisfies 

(id t + A D )vi = [A D , Xi}u, Vi\t= = Xi^(—^-) u o- (6.9) 

We introduce the Dirichlet Laplacian A^ acting on f2j := M 3 \ 0j and continue to denote 
At> the Dirichlet Laplacian outside the obstacle 0. Writing Duhamel's formula we get 

Vi(t) = e^xM^^o ~ i f e itA H(^)e- it ^ [A D , Xi}«r))dr. (6.10) 
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The solution Vi of (16.91) satisfies 



ll^ll^m £ «^-w«^l^ w + (6 ' u; 



* -A 



From [SI Thm.4.2] we know that the operator 

4-/<:= f T ^( Z ^)e- iTAD ^(r)dr (6.12) 



is bounded from L^H D 2 {Q) to Hp e (Q). Notice that in (16.121) we can introduce a cutoff 
function x% £ Co° equal to 1 on the support of Xi without changing the integral modulo 
smoothing terms. We need the next lemma: 

Lemma 6.2. Let ip e Qj°(R+) be a smooth function such that ip(r) = 1 if ip(r) = 1 and 
such that J2k>o^(^ 2kr ) = ^* e ^ e e Q U(J d to 1 on i/ie support of x%- Then for 

f e L'iSl) 

^(^UM^KXif) = Xi^G^){Xif) + Oj^njCA- 00 )*/. (6.13) 

We postpone the proof of Lemma 16.21 for the end of this section. 
End of the proof of Theorem \1.3l We introduce the operator 



■ ~ —A* 



which is continuous from L 2 (f2j) to L^H^iQ,^ with the norm bounded by A -Q / 4 . Indeed, 
since Xi,\ vanishes outside a small neighborhood of 0, we can apply Theorem 1 1.1 1 in M 3 \Gj. 
If we take fa = [Ajj, xi\ u -i then in view of Lemma [6.21 the last term in the right hand side 
of dnHD writes A i)X A*fi + L 2 {a) (\-°°). 

If Xi £ Co° ec i ua ls 1 on the support of Xi we can estimate 

n» e " A ^(^)«»n i|H | (n) < iw^)%ik< n) , 

where in the last inequality we applied [SI Thm.1.7]. Hence (16.111) becomes 

\\Xi,\ u \\ L2 H \ {Q) ~ A~7|| Xi A ^(^^) no || L 2 (n) + A-f +2e || XiiA ^(^^)uo|| L 2(n)- (6.15) 
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Set xx '■= J2i=iXi,\- Since {xi,x} have disjoint supports, (16.151) remains valid for Xi,x 
replaced by xx- We have thus obtained a smoothing effect with a gain a/4 — 2e and by 
interpolation with the energy estimate we have 

WxxuW^mn) < A" (1 - f+f) ||^(-^)u || LW (6.16) 
Away from <90 we can use the same arguments as in Section [5] and find 

11(1 - Xa)«||j4l.(0) < \ a \\*(-^-)u \\ LH n)- (6.17) 

Let p = 2 + e (for e > sufficiently small) and take ^(l — f + f) = a. 
Proof, (of Lemma \6.ty) : We fix A = 2 k °, k > 1. We write 

X^(-2" 2fc0 A D ) = ^ ^(-2- 2fc A J D )x^(-2- 2fc0 A D ) = (6.18) 

fc>0 

= ^(-2- 2fc0 A J D )x^(-2" 2fc0 A D ) + ]T ^(-2^ fc «(2^^ fc °)A^))x^(-2^ fc «A D ). 

Let B C M 3 be a ball of sufficiently large radius such that U^ 1 suppx« CC 5 and Bi C M 3 be 
balls such that suppxi C Bi. Suppose that Xt are suitably chosen such that the distances 
between any two such balls Bi be strictly positive (this is always possible, eventually 
shrinking the supports of Xi)- If we denote A^, resp. A l D the Laplace operators in the 
bounded domains B fl Q and resp. Bi fl fl, we notice that 

A D ( Xi f) = A D ( Xi f) = AU^f) = £$,(*</), v/ e L\fi). (6.19) 

Let 4> n , n > 1, resp. 0^, m > 1 denote an orthonormal system of eigenfunctions of A^, 
resp. A l D associated to the eigenvalues A 2 , resp. A^, i.e. such that — A D (p n = A 2 n , 
— A l D (j) l m = \^(p l m for n, m > 1. Consequently if /j G L 2 (f2) is supported in 5j fl fl we have 

= £ > ^(2- 2fe °A 2 )0 n , 

n>l 

and if g { G L 2 (f2) is supported in Bi (1 fl and we set A k = 2~ 2( - k ~ k °\ k > than 
^(-2- 2fc M fc A 4 D )^ = J2< 9i,tin > ^(2- 2fc M fe A^)^, 

m>l 

If the support of ^ is sufficiently large we show that the contribution of the sum in the 
second line of (16 . 1 8[) for k ^ ko is L 2^(X~°°). We consider separately the cases k > k , 
resp. k < k . Let first k > k : we distinguish two case, according to 2 ek ° < 2 h < 2 k ° for 
some e > 1/4 or 2 k < 2 k °' A . 
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• Let 2 k < 2 fc °/ 4 and set A k = 2 2( - k °~ k \ then A k > 2 k ° = A ~ A n . 

^(-2^ 2fe M fc A J D )x^(-2- 2fc0 A I) ) Xi / = (6.20) 

£ ^(2- 2k °A k \^(2- 2k °\ 2 n ) < Xi fAn >< XitnAln > C 
n,m>l 

< Xi0n, 4> l m >= T7 < -XiA D n , 0^ >= — < -Au(Xi^n) + [ A D, X#n, 0m > • 
K K 

(6.21) 

Since — A D (xi<fin) = — A^(xj0 n ) is self-adjoint, the first term in the last line writes 

1 ~ X i2 

and hence its contribution in the sum in (16.201) is 

^(2- 2k °A k \^(2~ 2k °\ 2 n ) < x^, <j? m > . (6.22) 

Since ip, ip are compactly supported away from 0, the only nontrivial contribution in 
the sum ( 16.201) comes from indices n, m such that A 2 ~ A k X£ ~ 2 2k ° and from ( 16.22ft 
this will be O^^^A^ 1 ) = Ol2^(1/ \ n ) which follows from the assumption 2 ek ° < 2 k 
for some e > 1/4. We estimate the last term in the right side of (I6.2ip 

< [A D ,x#„,C >= O^X- 1 ) = L2(n) (2- fc °), 

since on the support of ift, A n ~ 2 k °. By iterating these arguments M > 1 times we 
deduce that the contribution of the sum (16 . 18[) is Oi l 2^ j {2~ Mk °) for every M > 1. 

• Let now 2 efc ° < 2 k < 2 fc °, e > 1/4: in this case a simple integration by part is useless 
since the "error" is a multiple of the number of integrations. If the support of Xi is 
sufficiently small then (I6.19P holds and we have 

< Xi<frn, <f>ln >= J {Xi ~ l)0n</CA + 5 n=m , (6.23) 

where 5 n=m is the Dirac distribution. In the sum (I6.20p we see that the contribution 
from n = m is zero, since the support of ip(A k .) and are disjoint. For first 
term in the right hand side of (16.231) we use an argument of N.Burq, P.Gerard and 
N.Tzvetkov [9j Lemma 2.6]. Let k e 5(IR) be a rapidly decreasing function such that 
k(0) = 1. From a result of Sogge [T9], Clip. 5.1] we can write, on the support of \ 
where ( 16. 19ft holds 

«( <f-h D - X)f(x) = \ 1/2 J e lA ^a(x, y, X)f(y)dy + R x f(x), 
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where a(x,y,X) G is an asymptotic assumption in 1/A and —<p(x,y) is the 
geodesic distance between x and y, and where 

Vp, seN: H^a/II^CsuppxO < ^A- p ||/|| L2 . 



We use (16.191) and «;(v /= Ab - A n )0 n = n , k( v / -A' d - A^)0^ = 0^ in order to 
write < (xi — l)4> n , <fi l m > as an integral (modulo a remaining term, small) 



e 

x,y,z 



lXln/2Hx ' y ' z) a(x, y, X n )a(x, z, X i m )<p n {y)<^ m {z)dydzdx, 



with $(x, y, z) = ip(x, y) + a/ A^/ X^tp(x, z). Since |V z y?| is uniformly bounded from 
below and bounded from above together with all its derivatives, we obtain that 
there exist c > 0, Cp > such that |V X $| > c and <9^$ < Cp. It remains to 
perform integrations by parts in the x variable as many times as we want, each such 
integration providing a gain of a power Xn 1 ^ 2 - 



For k > k and 2 k ° < 2 fc / 4 , we write 

< -Xi(f>n, & D <P m >= T. 
A m A 



< Mm C >= T72 < -XiK ^ D 4L >=u2< -A D (xi0n) + [Ad, Xi]0n, 4L >, (6.24) 



in which case, using the spectral localizations ip, ip, we gain a factor from the first term 
in (I6.24p and a factor ^(A^) -1 from the second one; iterating the argument as many times 
as we want we obtain contribution L 2^(Al f ). In the last case 2 eh < 2 k ° < 2 fc , e > 1/4 
we use again the arguments of [HI Lemma 2.6]. The proof is complete. □ 



7 Appendix 

7.1 Hankel functions 

We will start by recalling some properties of Hankel functions that will be useful in deter- 
mining the behavior of the solution to (12. 2p (see [I]). The Hankel function of order v is 
defined by 

/+oo — itr 
e zsinht-ut dt ( 7J ) 
-oo 

For all values of u, {H u (z), H u (z)} form a pair of linearly independent solutions of the 
Bessel's equation z 2 y" + zy' + [z 2 — fi 2 )y = 0. 

Proposition 7.1. We have the following asymptotic expansions (see |2l Chp.9]): 
1. If v is fixed, bounded and z = rX ^> u > ^, - < 1 then 

H v {z) ~ ,/A e ^-T-f)(i + o^ 1 )), H u (z) ~ ,/A e -(Ar-^-f) (1 + 0(A -i )); 

V HZ V 7TZ 

(7.2) 
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2. If z = r\ > v >> 1 and ^ G [eo, 1 — ei] for some small, fixed eo > 0, e\ > 0, then 
writing - = —^-5 we have 

a v cos p 



H v {z)=H v {^-) ~ J—L^ e i ^f , -^(l + 0(y- 1 )l (7-3) 
cosp V ttz/ tan p 



H v {z) = H v {-^-s) ~ ,/_-— (1 + 0(z/ -i )); (74) 
cosp V nvtanp 

3. If z, v ^> 1 are nearly equal we have the following formulas 
(a) If z — v = tv x I z with fixed r, bounded, then 

J v {y + rv\) ~ ij^(-^2r)(l + O^ 1 )), (7.5) 



+ ri/*) ~ _^(5z(-^r)(l + O^" 1 )), (7.6) 



where for \r\ large and £ = |ra we /lave 

Az(-r) * ^TsinK+jKl+Ot 1 )), Si(-r) ^cosK+^l+Off 1 )); 

7T2T4 1 7T 2 T 4 

(7.7) 

If z = v, then 

2 1 / 3 

W "3^f(273T (1 + 0( )} ' (7 - 8) 
2 1 / 3 

nM ""3^f(273T (1 + 0(Z/ " )); (7 - 9) 
(c) If \v — r\\ < Co I -r A | then if vz = r\ we have for z < 1 (resp. z > 1) 

^)M^)-(^ + ^^o (1/ ^)), „ 10 ) 



where the function ( is defined by 

2 C 3/ 2 = f ^lEZ dt = ln[(1 + ^Y—2)/ z ] - VT^, z < 1, (7.12) 
3 J z t 



|(-£) 3 / 2 = / — l dt = v^T7-arccos(l/ 2 ), 2 > 1, (7.13) 
3 Jit 
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and where for r large and £ = |t 3 ^ 2 u>e /icrae 

Ai(r) = — Ue-*(1 + Or 1 )), 5t(r) = -^e«(l + Or 1 ))- (7.14) 

27T2T4 27T2T4 

Taking r = v 2/3 (, we compute Ai(u 2/3 () using (17.141) . (17.71) with £ = |( 3/2 ia 

Here J u and Y v are the Bessel functions of the first kind and H u (z) = J u (z) + iY u (z). 

Remark 7.2. We remark that ( defined in (17. 12ft . (17.1 3D is analytic in z, even at z = 1 and 
d(/dz < there; also, at z = 1, C = and (1 - z 2 )"^ = 2~ 2 / 3 . The formulas (ITTTOl) . fTTTTT]) 
are among the deepest and most important results in the theory of Bessel functions. These 
results do not appear in the treatise of Watson, having been established after Watson's 
second edition was published. See Olver [17] and [TJ. In the paper we use a simpler form 
of these asymptotic expansions for which we give an idea of the proof inspired from 

Proposition 7.3. For z > 1 close to 1 and v ^> 1 large enough we have 

V2 



J v {vz) ~ 

V 



%(-& i ) l,4 M-^C\ (7.15) 



Y v {uz) * -^(-^)^Bi(-u^(), (7.16) 

uniformly in z, where |C 3 ^ 2 = V z 2 — 1 — arccos(l/z) (see fl 7 ^). 
Proof. With a suitable choice of contour we have 

J v {uz) = — I ' e iw W' z) dt, (f>(t,z) = zsiat-t. (7.17) 
2vr J 

For v ^> 1 and z > 1 the saddle points are real, t = arccos(l/z) and the critical value is 
(p(t(z), z) = \/ z 2 — 1 — arccos(l/z). Now if ((z) is defined in terms of the exponent in the 
Debye expansion, it is analytic in z for z near 1 and we have 

Ai(-v 2 / 3 () = ±-J < '^^-^Ms = ^~J e^-X+Wdt (7.18) 

with critical points t 2 = £, thus we get (17.151) applying the stationary phase. We obtain 
the Debye approximations 

/ 2 Nl/2 

J v (vz 



Vz^l 



cosfz/V^ 2 — 1 — z/arccos(l/z) — 7r/4] (7-19) 



by replacing the Airy function by its asymptotic expansion, thus the proper condition for 
its validity is v 2 / 3 ^ ^> 1. For small u 2 ^ 3 ( we are in the regime v{z — 1) = rz/ 1 / 3 for which 
we have the estimations (17.51) . (17.61) . An extension of this result giving (17.101) . (17. lip uses 
a result of Chester, Friedman and Ursell who showed that a similar reduction is possible 
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whenever two saddle points coalesce: if d t (j)(t, 1) = and df t <f)(t, 1) = but df tt (fi(t, 1) ^ 0, 
then an integral of the form (I7.17P has a uniform asymptotic expansion in terms of the 
Airy function and its derivative. Their method was to make a change of variables so that 

<J ) (t J z) = (T-T 3 /3, (7.20) 

which holds exactly and uniformly; it is not merely an approximation for z near 1. □ 
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